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Abstract 

We argue that Kantor and Kardar's assertion that their simulation results 
contradict our criterion for the localization of a softly constrained ideal poly- 
mer is incorrect. Our criterion is inapplicable to the model used in these 
simulations. We argue that our criterion for localization of as polymer with 
internal constraints strongly suggests that the models with soft and hard 
constraints belong to different universality classes. The size measure used by 
Kantor and Kardar is also not adequate for studying polymer localization. We 
also show that, in spite of these contrasts, combining the simulation results 

with our criterion for ideal chains seems to uncover some intriguing physics. 
PACS numbers: 82.35. +t, 87.15.By, 61.41.+t 



1 



Typeset using REVTgX 



In the accompanying comment, Kantor and Kardar (KK) |T| discuss simulations of an 
ideal polymer chain of N monomers with M randomly chosen distance constraints. These 
simulations use an analogy between ideal polymers and resistor networks [[| to calculate 
the mean-squared end-to-end distance (r 2 ) of the polymer. The simulations together with 
a lower bound for r 2 are used to suggest that the polymer remains extended unless the 
number of constraints exceeds some M c which is of order N. They concluded that this 
result is at variance with our recent theory in which we argue that M c ~ N/ log N for an 
entirely different model (BT) ||. Here we show that this assertion is too strong, and that 
their results for ideal polymers may shed an interesting light on ours. 

The most obvious contrast between our study and that of KK is that the two studies 
investigate different models. In the random resistor analogy a short corresponds to a Dirac 
delta function (a hard constraint). In our Letter we explicitly excluded such delta function 
constraints; see footnote 9 of BT. Although a polymer with soft constraints will be more 
expanded than one with hard constraints, soft constraints add a new length scale (inducing 
longer range interactions) to the system which is present in our localization condition (see 
equations 1 and 10 of BT). The size of a polymer with soft constraints depends on both 
the number of constraints and the size of the constraints. Thus, the two models may 
be in different universality classes, making the objections of KK circumspect. The above 
notwithstanding, expanding our localization condition (eq. 10 of BT) for large = {b/c) 2 
and small 5 = 1 — /i gives the localization condition 

log + log L — 1 

The hard constraint case corresponds to the limit — > oo with ^ > £ > 1, yielding a 
maximum value of 5 — > 0, self-consistent with our expansion and consistent (!) with the 
arguments of KK that the critical value of M ~ N in this limit. 

The purpose of our Letter was to determine whether a polymer with soft internal con- 
straints is better approximated by an expanded chain or by a collapsed chain. The nature 
of the transition between these two regimes could not be determined by our methods . A 



smooth, continuous transition, as found in the simulations of KK is not inconsistent with 
our results. The question then becomes, what is the size of the polymer at this transition? 
If we assume that R G oc N/M holds for the radius of gyration (Rg) of our model then 
at the transition between the two regimes, R G ~ logiV. This scaling is a reasonable divid- 
ing line between the fully expanded (R G ~ N) and fully collapsed (R G ~ 0(1)) regimes. 
Furthermore, this dividing line is not arbitrary. It is consistent with a completely differ- 
ent criterion for polymer localization, as we now show. For a general system of N (ideal) 
monomers connected by harmonic bonds, 



where the Ui are the non-zero eigenvalues of the connectivity matrix, w min is the lowest 
non-zero eigenvalue, u) max is the largest eigenvalue, and n{uS) is the spectral density H0. 
The eigenvalue uj min is related to large scale properties of the network, whereas u max is 
related to local, non-universal properties of the network ||. If we take the network to be 
fractal, then n(uj) ~ uj^ 3 " 2 ^ 2 , where d s represents the effective spectral dimension of the 
network ||. An unconstrained ideal polymer has d s = 1 and d s increases with the number 
of constraints. For d s > 2, R G depends on u max , and therefore on the local properties of 
the network. For d s < 2, u min ~ N~ 2 / d % so R 2 G ~ uj { ^~ 2)/2 ~ N 1 - 2 -^^ §. The transition 
between these two regimes occurs when d s — 2 and therefore R G ~ logA^, which is precisely 
the dividing line obtained from using our localization criterion with R G oc N/M. Although 
this derivation approximates the constrained chain as a homogeneous fractal, it definitely 
shows that, if the results of KK and of Solf and Vilgis 0] are correct for our model, then the 
localization criterion is related to the behavior of general Gaussian networks in an intriguing 
and unexpected way. 
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